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ABSTRACT 

This work outlines a procedure for simulating the 
flow field within multistage turbomachinery which 
includes the effects of unsteadiness, compress i bi 1 i ty , 
and viscosity. The associated modeling equations are 
the average passage equation system which governs the 
time-averaged flow field within a typical passage of a 
blade row embedded within a multistage configuration. 
The results from a simulation of a low aspect ratio 
stage and one-half turbine will be presented and com- 
pared with experimental measurements. It will be 
shown that the secondary flow field generated by the 
rotor causes the aerodynamic performance of the down- 
stream vane to be significantly different from that of 
an Isolated blade row. 

INTRODUCTION 

A goal of computational fluid dynamics for turbo- 
machinery is the prediction of performance parameters 
and the flow processes which set their values. Achiev- 
ing this goal for multistage devices is made difficult 
by the wide range of length and time scales in the 
associated flow fields. Currently the procedure used 
in design and off design analysis is based on a quasi- 
three-dimensional flow model whose origins can be 
traced back to the late forties and early fifties 
e.g., <Wu, 1952; Smith, 1966). This model requires 
calculations to be executed on two orthogonal surfaces 
within a blade row passage of a multistage configura- 
tion. One of these surfaces is an axisymmetric sur- 
face of revolution whose intersection with a blade row 
defines a cascade. The flow field relative to this 
cascade is assumed to be steady in time. In practice, 
a finite number of such surfaces are chosen to define* 
a series of cascade flows from hub to shroud. The 
other surface represents a meridional through-flow sur- 
face. The flow associated with this surface is an axi- 
symmetric representation of the flow field within the 
machine. This flow field is also assumed to be 
steady. The flow fields on both surfaces are coupled 
and are solved iteratively. The effects of unsteadi- 
ness, turbulence, and endwall secondary flows are 
Introduced through empirical correlations. 


Although proven to be very useful, this flow 
model has its limitations. Among these is off-design 
performance analysis, and the ability to analyze uncon- 
ventional machinery where extrapolation of the underly- 
ing empirical database is required. Other problems 
arise whenever there are large local variations in the 
radial velocity component within a blade passage. 

Such variations can be brought about by in-passage 
shock waves, separated boundary layers, and endwall 
secondary flows. It is generally agreed upon that a 
way of overcoming these shortcomings is the develop- 
ment of a true three-dimensional flow model. 

Two three-dimensional flow models have been pro- 
posed for the simulation and analysis of multiple 
blade row flows. The first (Denton, 1979; Adamczyk 
1984; N1 , 1987) referred to as the average passage 
flow model in Adamczyk (1986) simulates the time- 
averaged flow field within a typical passage of the 
blade row. The second simulates the unsteady determin- 
istic flow field within the machine. Although a 
number of unsteady simulations of single stage turbine 
configurations and counterrotating propellers have 
been reported (Rai, 1987; Whitfield et al., 1987), 
executing an unsteady simulation of a multistage con- 
figuration of practical interest is far beyond the 
capabilities of todays advanced supercomputers. Fur- 
thermore, it is by no means obvious that performance 
prediction requires such a high degree of flow resolu- 
tion. However, because unsteady simulations of an 
idealized configuration may prove to be a useful means 
of investigating the closure issue associated with 
time-averaged flow models, this activity should be pur- 
sued. In this work it will be shown that the simula- 
tion of the time-averaged flow field within multistage 
machinery is within the capabilities of todays advanced 
computers and that the average passage flow model 
gives more insight into the flow phenomena that con- 
trol the performance of multistage machinery than 
today's quasi three-dimensional flow models. 

The objective of this paper is to outline a proce- 
dure for simulating the time averaged flow field within 
a typical passage of a blade row within a multistage 
machine. This model includes the effects of viscosity 
and compressibility, and the influence of neighboring 
blade rows. The mathematical formulation upon which 


ORIGINAL PAGE IS 

OF POOR QUALITY 


this model is based has been outlined in Adamczyk 
(1984). The algorithm used to solve the inviscid form 
of the governing equations is reported in Celestina, 
(1986); Adamczyk et al. (1986). The current work out- 
lines a numerical solution procedure for the viscous 
form of these equations and an acceleration technique 
to enhance convergence. In addition, a comparison 
will be made between experimental data recorded during 
tests of a one and one-half stage, large scale, low 
speed, axial flow research turbine and simulation pre- 
diction. The underlying unsteady flow physics which 
appears to control the performance of the second vane 
of this machine will also be discussed. 

GOVERNING EQUATION 

A complete derivation of the three-dimensional 
average-passage equation system is presented in 
Adamczyk (1984). These equations were derived by 
filtering the Navier Stokes equation in both space and 
time to remove all information except that associated 
with the time average flow field within a typical 
passage of a blade row of a multi-stage configuration. 
With respect to this blade row, the integral form of 
these equations can be written 


J (Xq)dV + L(Xq) = J XS dV + 


Xk dV 


L v <Xq> 


(1 ) 


The vector q contains variables density, axial and 
radial momenta, angular momenta, and total internal 
energy. X is the neighboring blade row blockage fac- 
tor and ranges between zero and unity, unity being the 
value associated with zero blade thickness. This 
parameter explicitly introduces the effect of the 
neighboring blade row blade thickness. The operator 
L(Xq) balances the mass, axial and radial momenta, 
angular momentum, and energy through a control volume. 

Xk dV is a source term due to the cylindrical 


coordinate system and 


XSdV contains the body 


forces, energy sources, momenta, and energy temporal 
and spatial mixing correlations associated with the 
neighboring blade rows. A procedure for estimating 
S has been outlined in Adamczyk et al . (1986) and is 
extended here to include the effects of viscosity. 

The operator L v (Xq) contains the viscous and heat 
transfer terms. The vector q and the operators L 
and are defined as 
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L = | CXF dA + XG dA + XH dA n ] (3) 

dA z rQ 


In the above equations p represents the density, v 
the absolute velocity vector, p the pressure, and T 
the temperature. The differential dV is the volume 
of the control volume and dA z , dA r , dA 0 are the dif- 
ferential areas of its sides. From the equation of 
state, the total internal energy is related to pres- 
sure through the equation 


and 


L » - U 1 


L = [XF k dA, + XG dA + XH dA-] (4) 
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and the total enthalpy, H, is related to p and e 0 
by 


where 


— 2 T 

F = Cpv 2 .pv 2 + P,pv 2 v r r p v 2 v Q ,pHv z ] 

- 2 T 
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( 6 ) 


H = e„ 


( 22 ) 


Sutherland's law is used to determine the molecular 
viscosity coefficient, and Stokes' hypothesis gives 
X v = -2/3 pg,. Turbulence is accounted for by adding a 
turbulent viscosity pt to the molecular viscosity pq. 
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+ Mt 


(23) 


In a similar manner, the molecular thermal conductiv- 
ity k is replaced by 


k = C 


4L 


(24) 


where Cp is the specific heat at constant pressure 
and P r | <i + P r 1 1 is the laminar and turbulent Prandtl 
number respectively. The two-layer algebraic model of 
Baldwin and Lomax (1978) is used to model pf 

All lengths in the above equations are nondimi- 
sionalized by a reference length normally taken as the 
largest blade row diameter. The velocity components 
are nondims ional ized by a reference speed of sound, 
a ref^ Y wh ere y Is the ratio of specific heats. 
Pressure and density are nondi mensional i zed by their 
respective reference values. 

For rotating flows, the absolute (fixed) refer- 
ence frame is transformed to the relative (rotating) 
frame by the transformation 


Absolute = ^relative + 

where Q is the rotational wheel speed. 
Eq. (24) into Eq. (1) transforms L and 


(25) 


Introducing 
L u to 


- u 


XF dA + 


XG v dA f 


\(H - rQq)dA 


(26) 


and 


L v = 


JdA^v 


dA. 


XG v dA r 


X( H 


rQq) dA e J (27) 


Discretization of the inviscid portion of Eq. (1) 
including estimates of the surface area and volume is 
presented in Celestina et al . (1986). The viscous and 
heat transfer portion of Eq. (1) is discretized by 
evaluating the shear stress and the heat flux at the 
center of each face. The shear stresses and heat flux 
are estimated using central differences of the velocity 
and the temperature field. 

ARTIFICIAL DISSIPATION 

To suppress odd-even point decoupling of the solu- 
tion to the discretized equations, dissipative terms 
are added to the equations. The operator L v is 
replaced in the discretized form of Eq. (1) with D(Xq) 


D(Xq) » Dj(q) + L v (Xq) 


(28) 


where Dj ( q ) is the added artificial dissipation opera- 
tor needed to prevent decoupling of the solution within 
inviscid regions of the flow. The operator Dj is 
patterned after the model developed by Jameson et al . 
(1981) and is composed of three spatial operators. 


Dj(q) - (0 2 + D r + D 0 )q 


(29) 


which can be evaluated separately. The dissipation 
in, for example, the axial direction (and similarly 
for the others) is expressed as follows: 


where 
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second difference operators, while 13 is the maximum 
eigenvalue of the Jacobian matrix formed from F. 

The coefficient vj j k is defined as 


(4) are constants set at 1/8 and 1/512 respec- 
The symbols A and A^ denote the first and 


i.j.k 


P i^l,.j,k ZP i.i,k 


P. 


L-.I.j ,k 


P i + l,j.k + 2P i,j,k + P i-!,j,k 


(34) 


and is used primarily to prevent oscillations near 
stagnation points and shocks. 

Vatsa and Wedan (1988) scaled the Jameson artifi- 
cial dissipation operator by a function of the local 
Mach number to reduce its effect within viscous 
regions of the flow. The present work uses the local 
meridional Mach number M normalized by a reference 
upstream meridional Mach number M to accomplish this 
task. To prevent this function from increasing the 
level of dissipation in the inviscid flow regions, the 
maximum value of this function is taken as unity/ To 
reduce. the level of artificial viscosity resulting 
from highly stretched mesh cells, we also adopted what 
is referred to by Vatsa and Nedan as individual eigen- 
value scaling of the artificial dissipation operator. 

SOLUTION PROCEDURE 

The discretized form of the equations associated 
with the averaged passage flow model are solved using 
a dimensional sequencing algorithm. The motivation 
for the present algorithm came from observing the 
evolution of the error history associated with the 
algorithm reported in Celestina (1986). For many 
cases the magnitude of the error associated with the 
axisymmetric component of a variable was a significant 
fraction of the magnitude of the error associated with 
the variable itself. It thus seems reasonable to 
expect a solution algorithm which explicitly reduced 
the error of the axisymmetric component of the flow 
field would enhance the convergence of the three- 
dimensional field. With the body forces assumed 
known, (i.e., S assumed known) the present algorithm 
iterates between the three-dimensional flow equations 
and the corresponding through-flow equations to enhance 
the rate of convergence to the three-dimensional time 
asymptotic flow problem. The construction of this 
algorithm is as follows. First, the through-flow equa- 
tions compatible with Eq. (1) are derived by summing 
Eq. (1), as modified according to the discussion in 
the preceding sections, over the tangential index k. 
This is accomplished by premultiplying Eq. (1) by the 
operator 

Sfm i L (35) 

k*l 
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where K Is the number of control volumes spanning 


the pitch. The result may be written as 

Xq dV ♦ L(q) - D(q) - J Xk dV + S dv + Xk dv 

- J X k dv + L(ifq) - L(q) - D(ifq) +i"D(q)| (36> 

where dv, X, q are defined as (37) 

ifX dV - X dv (38) 

ifX dv - Xq dv (39) 


while the operators if , L and D are 


if q - — X g -i V -~q (40) 

X dV 


L(q) - if L(q) (41) 


D(q) - ifD(q) (42) 

The terms which appear on the right hand side of 
Eq. (36) are treated as forcing functions and are esti- 
mated using the most recent value of the three- 
dimensional flow variables. Note that the expressions 
which appear within brackets vanish by construction In 
regions of the flow where q * q, and upon convergence 
of Eq. (36) (i.e. 8/9t Xq dv * 0)^q * q. The time 
asymptotic solution of Eq. (36) is thus identical to 
the axlsymmetrlc average of the time asymptotic solu- 
tion to Eq. (1). The steady state solution of Eq. (1) 
can thus be obtained by cycling between a time- 
advancing algorithm for Eq. (1) and a similar algo- 
rithm for Eq. (36) . 

The discretized form of Eq. (1) is advanced In 
time using the four-stage Runge Kutta algorithm of 
Jameson et al . (1981). Local time stepping (i.e., con- 
stant C.F.L. number) and residual averaging Is 
employed to enhance convergence. Upon completing of a 
fixed number of temporal relaxation cycles, the three- 
dimensional flow variables are used to evaluate the 
right hand side of Eq. (36). This equation Is then 
advanced In time using the same Integration procedure 
as that for the three-dimensional system. After a 
fixed number of time steps, the value of q Is updated 
according to the equation 

q - q + (q - <£ q) (43) 


where the q's within the brackets are those that were 
used to evaluate the right from hand side of Eq. (36). 

The three-dimensional residual error based on the 
updated value of q (I.e., Eq.( 43)) Is generally 
found to be largest within the blade passage region. 
Prior to Initiating the next three-dimensional Itera- 
tion cycle, the residual Is reduced by performing a 
fixed number of three-dimensional iteration cycles 
(between five and ten Iterations) over a local three- 
dimensional flow region which extends a modest dis- 
tance upstream and downstream of the blade passage 


leading and trailing edges. The Inlet and exit bounda- 
ries of this reduced flow domain are located in a 
region where the sum of the bracketed terms In Eq. (36) 
Is small. The Inlet and exit boundary conditions asso- 
ciated with the reduced flow domain are derived using 
the updated variables obtained from Eq. (43). Within 
this region of flow the values of q obtained from 
relaxing the reduced flow domain equations replace 
those obtained from Eq. (43). 

The outlined solution procedure is a multigrid 
algorithm of the form Introduced by Brandt (1982). It 
was specially constructed to reduce the axlsymmetrlc 
component of the error vector, which at the start of 
the solution procedure Is often large. It also recog- 
nizes the spatial relaxation of the three-dimensional 
flow field to an axlsymmetrlc field away from the 
blade row of Interest. This recognition reduces the 
computational work required to obtain a solution rela- 
tive to a more traditional multigrid strategy. 

The solution procedure outlined above has been 
Implemented In both a V and H cycle framework. In the 
V cycle strategy, one proceeds directly from the 
three-dimensional solver to the through-flow solver to 
the reduced flow solver before going back to the 
three-dimensional solver. In the W cycle, one pre- 
ceeds from the three-dimensional solver to the 
through-flow solver to the reduced flow solver and 
then cycles between the through-flow solver and the 
reduced flow solver before going back to the three- 
dimensional flow solver. The simulation to be reported 
was executed using the V cycle strategy. Experience 
with the W cycle Is limited and needs further develop- 
ment. However, a preliminary analysis showed the W 
cycle strategy to require less computational work than 
the V cycle to converge to a fixed tolerance level. 

When the three-dimensional flow field converges 
to a predetermined level, the body forces and energy 
sources required as Input to simulations of the remain- 
ing blade rows can be estimated using the procedure 
outlined In Ref. 9. The cycling of Information 
between the simulated blade rows of the multiblade row 
machine Is carried out until the tangential average of 
each simulated blade row flow field agrees with each 
other to a predetermined tolerance. 

BOUNDARY CONDITIONS 

All solid surfaces are modeled as rigid, nonslip, 
and Impermeable. The surfaces are also assumed to be 
adiabatic. These conditions imply that the velocity 
relative to a solid surface is zero and that the tem- 
perature gradient normal to the surface Is also zero. 

The pressure at a solid surface Is obtained from the 
normal momemtum equation evaluated at the surface. At 
the Inlet, either the mass flow or the total pressure 
is specified along with the total temperature and the 
radial and tangential velocity components. The one- 
dimensional Relnmann invariant C” Is extrapolated from 
the Interior to the boundary; with the specified flow 
variables. It defines the Incoming pressure, axial 
velocity component, and temperature. The shear 
stresses and heat flux at the inlet are also set to 
zero. At the exit radial equilibrium, with the pres- 
sure specified at the hub, Is used to establish the 
radial pressure distribution. The flow quantities p, 
pv z , pv r , pv e are extrapolated from the Interior. 

GRID GENERATION 

As discussed in Celestlna et al. (1986), the aver- 
aged passage equation system requires that a mesh be 
specified for each blade row of a multistage machine. 

In addition, the meshes must have a common meridional 
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mesh in order to eliminate the need for interpolating 
the body forces and correlations from grid to grid. 

To capture shear layers and stagnation points properly 
a fine mesh spacing is required In a direction normal 
to solid surfaces and in the blade leading and trail- 
ing edge regions. A mesh generator which is capable 
of generating these features is discussed in detail in 
Mulac (1986). 

RESULTS AND DISCUSSION 

The simulation executed was the Low-Speed 
Rotating Rig. at United Technologies Research Center. 
The Low-Speed Rotating Rig (LSRR) is a stage and a 
half turbine consisting of an inlet guide vane, a 
rotor, and a stator. The inlet guide vane contains 22 
blades and the rotor and stator both contain 28 
blades. The flow coefficient, <J> , Is 0.78 and the 
spacing between blades, B x , is 0.5. The LSRR grid 
contains 228 axial, 25 radial, and 41 ci rcumferenti al 
points. Each blade row contains 40 axial points dis- 
tributed along the chord with 26 axial points between 
each blade row, the inlet and exit. 

The results to be presented required eleven hours 
of Cray 2 C.P.U. time. They represent but a small 
fraction of the information obtained from the 
simulation. They are intended to illustrate the degree 
to which one can quantitatively predict performance 
parameters which are of interest to designers, and to 
reveal qualitative information identifying flow phenom- 
ena which may have an impact on performance. These 
results also reflect the current state of model devel- 
opment. The first series of results shows the pre- 
dicted pressure distribution on the surface of each 
blade row of the turbine as a function of axial chord 
length and percent of span height. The span locations 
measured from the hub are 1.3, 12.5, 50, 87.5 and 98.7 
percent, respectively. The experimental measurements 
taken at these locations are also shown. Experimental 
data was also available for 25 and 75 percent of span 
but was not utilized since it provided little addi- 
tional information relative to the current discus- 
sion. The results for the first vane are shown in 
Fig. 1. The predicted loading level is seen to be in 
good agreement with the measurements of Dring (1988). 
The predicted pressure surface pressure distribution 
is in excellent agreement with the experimental 
results. For the suction surface, the agreement 
between measurement and simulation is good for the 
region forward of the minimum pressure peak. Aft of 
the peak, the agreement between experiment and simula- 
tion deter iorates . This deterioration is believed to 
be related to viscous effects (i.e., turbulence and 
transition modeling) whose modeling could be improved. 
Some exploratory calculations suggest that the boundary 
layer aft of the suction surface minimum pressure is 
growing too rapidly and, as a result of the radial 
pressure gradient, is being transported towards the 
hub to an extent greater than that suggested by a flow 
visualization studies. Improvements in the agreement 
between simulation and experiment have been obtained 
by incorporating a simple transition model in which 
the flow remains laminar forward of the minimum pres- 
sure peak and Baldwin Lomax turbulence model as imple- 
mented by Dawes (1986). The current simulation assumes 
the flow to be fully turbulent from the leading edge 
of each blade. Work on this problem is continuing. 

The next figure shows the predicted and measured 
pressure distribution for the rotor. The predicted 
loading levels appear to be in good agreement with 
measurements, with the exception of the hub and tip 
region. The present simulation does not include a 
clearance region, which should account for some of the 


discrepancy in the tip region. The pressure distribu- 
tion along the pressure surface is once more in excel- 
lent agreement with the measurements. At the midspan 
and at 25 and 75 percent (not shown) of span the pre- 
dicted pressure distribution along the suction surface 
is In good agreement with the data. At 1.3 percent 
and 12.5 percent of span, the suction surface pressure 
coefficient Is lower than that measured. As a result 
the loading is lower over the forward portion of the 
rotor than what has been measured. Although the cause 
of this discrepancy is unknown at the present time, 
one could speculate that it may be due to an over- 
estimate of the magnitude and extent of the low momen- 
tum fluid exiting the first vane. 

The pressure distribution for the last vane is 
shown in Fig. 3. Once again the loading level is 
well predicted with the exception of location that is 
1.3 percent of span. The underpredicted suction sur- 
face pressure coefficient at 1.3 percent of span. The 
underpredicted suction surface pressure coefficient at 
1.3 and 12.5 percent of span suggests that the flow 
incidence to these sections is underestimated. There 
also appears to be a shift of the predicted pressure 
distribution relative to the measured distribution. 

This shift is believed to be caused by an over esti- 
mate of the loss generated by the first two blade 
rows. With the exception of this discrepancy, the 
pressure distribution on the pressure surface is in 
good agreement with measurements. Similarly, the pre- 
dicted suction surface distribution at midspan agrees 
well with the experimental distribution. The next fig- 
ure shows the predicted relative total pressure coeffi- 
cient forward and aft of the rotor as a function of 
span. The measured distribution reported in Sharma 
et al. (1988) is also shown. The magnitude of the pre- 
dicted coefficient for the inlet flow is higher than 
measured; however, the shape of the curve is consist- 
ent with the data. The magnitude of the predicted 
exit flow coefficient is also higher than measured. 

The Influence of the secondary vortices generated 
within the rotor passage on the exit flow coefficient 
is more evident in the experimental data than in the 
simulation result. The data in Sharma et al . (1938) 
suggests that the secondary vortices exit the rotor at 
approximately 30 and 70 percent of span. The local 
minimums in the measured exit flow distribution at 25 
and 85 percent of span are a consequence of the veloc- 
ity field induced by these vortices. The location of 
the tip vortex as suggested by the experimental data 
is significantly inboard of the location 90 percent 
span suggested by the simulation. The difference is 
believed to be due to the tip leakage flow which was 
not accounted for. It appears that this flow drives 
the tip secondary vortex inward towards the hub. 

The simulation results place the hub secondary 
flow vortex at 25 percent of span at the exit of the 
rotor, which is in good agreement with measurements. 
This vortex, however, appears to be more diffuse than 
measured, which may account for the lack of a local 
minimum in the exit flow relative total pressure coef- 
ficient at 25 percent of span. 

The deviation from the design Intent of the rotor 
relative exit flow angle is shown in Fig. 5. This 
angle is plotted as a function of span, with the solid 
curve representing the simulation and the open symbols 
the measured data from Sharma et al , 1988. The agree- 
ment between the two is reasonable with the exception 
of the tip region, where the result of neglecting the 
tip leakage flow is quite evident. The noticeable 
overturning of the flow near the endwalls and the sub- 
sequent underturning in the region of midspan caused 
by the secondary vortices is clearly seen in both the 
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data and the prediction. The flow physics which leads 
to this result appears to be well captured. 

Recently the authors became aware of data 
acquired under A.F.O.S.R. sponsorship by United 
Technologies Research Center which shows the agreement 
between simulation and measurement to be better than 
that indicated by Figs. 4 and 5. An evaluation of the 
present results in light of this data will have to 
await its publication. 

The next two figures are for the second vane. 

The first, Fig. 6, shows the predicted total pressure 
coefficient forward and aft of the second vane as a 
function of span. Also shown is the data from Sharma 
et al . (1988) . 

The predicted total pressure coefficients are lower 
than measured which, as previously noted, is the 
result of overestimating the loss produced by the 
first two blade rows. The difference between these 
two curves is a measure of the loss across the blade 
row. It too appears to be overestimated. With the 
exception of the flow region strongly influenced by 
the tip leakage flow from the rotor, the indicated 
trends agree with the experiment. Judging by the 
shape of the incoming total pressure field (either 
from the data or the prediction) it does not appear 
that the secondary vortices generated by rotor would 
have a significant effect on the performance of the 
second vane. The next figure, which shows the flow 
angle entering and leaving the second vane relative to 
the design intent, leads to a far different conclu- 
sion. The large variation in the inlet flow angle 
between 10 and 40 percent of span seen in both the 
experimental data and prediction clearly is caused by 
the rotor hub secondary vortex. The time-averaged sig- 
nature of this vortical structure is a region of shear 
whose vorticity vector is nearly aligned with the 
incoming time-averaged velocity field. The tip second- 
ary vortex appears to generate a similar structure; 
however, the measurements suggests it occurs in board 
of that predicted. This entering shear flow causes 
the second vane to behave aerodynami cal ly differently 
than it would in isolation. 

The results shown in Fig. 7 indicate that the sim- 
ulation predictions are in reasonable agreement with 
data. The simulation captures the underturning of the 
flow near the endwalls, and in the midspan region, the 
tendency for the flow to overturn. This behavior is 
part due to the time-averaged shear flow produced by 
the rotor secondary vortices. 

The last figure shows the limiting streamline pat- 
tern associated with the time-averaged flow adjacent 
to the second vane pressure surface. Also shown is 
the corresponding flow visualization result of Lang- 
ston (Langston, 1988). These results are presented to 
show that the time-averaged effect of the unsteady vor- 
ticity field entering a blade row has a major influ- 
ence on the flow field within the blade row. Near 
midspan, the incoming shear flow causes a contraction 
of the streamlines, suggesting the existence of a line 
of flow separation. Near both endwalls, the time- 
averaged vorticity field entering the second vane 
interacts with the secondary vortices generated within 
the vane to generate two lines of flow attachment. 

This flow pattern is far more complicated than that 
for an isolated blade row and the ability of the 
present model to capture its structure is very encour- 
aging. Finally, in closing it is suggested that the 
enhanced heat transfer observed at midspan of the suc- 
tion surface reported in Sharma et al. (1988) is 
caused by the secondary vortices generated within the 
rotor . 


SUMMARY AND CONCLUSION 

Given the early state of the average passage 
model development, the results presented in this 
report are very encouraging. The amount of empirical 
information used in the stage and one-half turbine sim- 
ulation is considerably less than that required to 
achieve comparable results using today's quasi three- 
dimension flow models. The average passage flow model 
also appears to be able to capture the physical fea- 
tures of the secondary flow field generated within a 
multistage turbine configuration. It appears that an 
accurately modeled of the time-averaged vorticity 
field produced by the unsteady secondary vorticity 
exiting a blade row is required to establish the per- 
formance of the downstream blade row. The outlined 
closure procedure, unlike some others which have been 
suggested, insures that the unsteady vorticity field 
exiting a blade row is consistent with the exiting 
time-averaged vorticity field. Hence, no spurious vor- 
ticity is produced as a result of coupling one blade 
row to another. 

The simulation of the stage and one-half turbine 
has shown the complex nature of the endwall flow in a 
low aspect ratio turbine. In the rotor, the secondary 
vortices generated within the endwalls affect the flow 
at midspan. They cause the endwall fluid to deposited 
on the suction surface. This transport of low momen- 
tum fluid leads to significant spanwise mixing across 
the axisymmetric stream surface. It is doubtful that 
an endwall boundary layer model could be used to pre- 
dict this phenomenon. 

There are numerous research activities that need 
to be pursued to further the development of the aver- 
age passage flow model. Other configurations need to 
be simulated, including multistage compressors and 
high- speed machinery. Grid generators need to be 
developed which are compatible with the average pas- 
sage model that can cluster grid points in regions of 
high flow gradients that occur near and away from solid 
surfaces. Algorithm improvements need to be made to 
speed convergence. The sensitivity of design parame- 
ters to turbulence modeling for the average passage 
model must be established. Finally, experimental and 
analytical work in support of closure modeling must 
be pursued to establish the models for generalized 
Reynolds stresses and the associated energy 
correlations. 
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